Quantum Zeno-like effect due to competing decoherence mechanisms 
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We propose a selfconsistent quantum mechanical approach to study the dynamics of a two-level 
system subject to random time evolution. This randomness gives rise to competing effects between 
dissipative and non-dissipative decoherence with a consequent slow down of the atomic decay rate. 
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I. INTRODUCTION 
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Time in Quantum Mechanics plays a rather ambigu- 
ous role. Usually, it enters as a continuous parameter, 
but the principle of relativity suggests a parallelism be- 
tween position- momentum and time-energy. Thus, there 
was the attempt to quantize the time |l|]. Otherwise, 
it has been considered as a discrete variable (2||^. This 
approach implies a modification of the Schrodinger equa- 
tion providing also an explanation of the nonappearance 
of macroscopically distinguishable states in terms of non- 
dissipative decoherence. 

More recently, a generalization of the Liouville-Von 
Neumann equation [Q was developed without any specific 
statistical assumption . It results an expression for the 
density operator which, a posteriori, can be interpreted 
as if the evolution time is not a fixed parameter, but a 
stochastic variable whose distribution is a F-distribution. 
It also provides a model for quantum measurement. Fur- 
thermore, the characteristic time appearing in such ap- 
proach could be related to the time-energy uncertainty 
relation |^,||. This approach is also useful to explain 
decoherence whenever the environmental degrees of free- 
dom responsible for decoherence are not easily recogniz- 
able 0. 

Here, we shall adopt this approach for a dissipative 
two-level system. Then, we shall show how decoherence 
effects compete in this framework in different regimes. 
In particular, we shall give an example of how statisti- 
cal quantum theory can describe a Zeno-like effect, i.e. 
frozen decay, which is generally thought to arise from the 
state reduction caused by the measurements M . 



II. STANDARD APPROACH 

Usually, the spontaneous decay of a two-level atom is 
described by means of a master equation derived by as- 
suming the system interacting with an environment ^ . 

Consider a atom, with two relevant levels {\g) , |e)} and 
lowering operator a = |.g)(e|. Let u be the energy dif- 
ference between the two levels {Ti = 1), and let 7 be the 
decay rate. Then, the master equation we are concern- 
ing is derived by considering an interaction of the type 
H oc (crFt + crtr) where F is a bath operator. Under the 
Born-Markov approximation one arrives to |10[| 



with the Liouvillian superoperator 

Cp — —iuj [oztP] + — {2apa* — a^ap — pa* a) , (2) 

where, we choose to define the inversion operator as cr^ = 
cr^cr — aa'^ , and the quadrature operators as — <y ^ cr'^ 
and <Ty — ia — ia'^ . 

Dephasing processes, if necessary, are introduced in the 
same way. These processes do not change the population 
of the two-level atom but do cause a phase randomiza- 
tion of the atomic dipole. Then, they can be modeled 
by considering the inversion ct^ to be coupled to environ- 
ment by oc crz(F + F^). The master equation (|^) now 
becomes 



p = -iuj [0-2, p] 
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where k represents the phase decaying rate. 

From the master Eq. (||) it is easy to derive the follow- 
ing dynamical equations 



Tr {pa^} = -7 Tr {pa^} - 7 , 



TrW 



(4) 
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The solutions read 



Tr{p(<)(7j =Tr{p(0)aJexp(-7i) 
+ [cxp {--ft) - 1] , 

TT{p{t)a} = Tr{p(0)cr}exp \-iujt~ 



4r + K]t 



(6) 
(7) 



We may see that the equation of motion for the inver- 
sion is unchanged with respect to the dissipative case, 
but the polarization decay rate is increased above the 
spontaneous emission result. 

Nevertheless, decoherence is not always necessarily due 
to the entanglement with an environment but it may be 
due, especially the non-dissipative one, to the fluctua- 
tions of some classical parameters or internal variable of 
the system. Or it might have an "intrinsic" character 
pd| . Hence, we shall present a more general approach to 
non-dissipative decoherence for a two-level system. 
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III. RANDOM TIME EVOLUTION 

Quantum mechanics is a statistical theory whose ele- 
ments are ensembles of quantum systems, or ensembles 
of measurements on the same quantum system. This led, 
long time ago, to the introduction of the density operator 
Q. Along this line, we cannot state a priori that time 
is uniquely determined within the ensemble. Rather, it 
would be more reasonable to give a statistical interpre- 
tation of the time variable too. Then, following Ref. 
the evolution of a system is averaged on a suitable proba- 
bility distribution P{t, t') where t' represents all possible 
times within the ensemble. Let p{0) be the initial state, 
then the evolved state would be 



dt' P{t,t') p{t') 



(8) 



where p{t') — exp{— i£t'}/9(0) is the solution of the 
Liouville-Von Neumann equation Q. 
One can write as well 



-p{t) = V{t)p{Q) , 
where the superoperator V is given by 



V(t) = / dt' P{t, t') e 



-iCt 



(9) 



(10) 



In Ref. 1^, the function P(t^t') has been determined 
to satisfy the following conditions: i) 'p{t) must be a 
density operator, i.e. it must be self-adjoint, positive- 
definite, and with unit-trace. This leads to the condition 
that P{t,t') must be non-negative and normalized, i.e. a 
probability density in t', so that Eq.(^ is a completely 
positive map; ii) V(t) satisfies the semigroup property 
V{ti +t2) = V(ti)V(i2), with ti,t2 > 0. These require- 
ments are satisfied by 
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and 



P{t,t') 



(1 + iCry/^ 
tWM 



(11) 



(12) 



where the parameter r naturally appears as a scaling 
time. Notice that the evolution superoperator ( [ll] ) only 
depends on t, and parametrically on r as in "non exten- 
sive" generalization of Liouville equation Indeed, t' 
comes out when a statistical interpretation is employed. 
Expression ( |l2| ) is the so-called F-distribution function, 
well known in line theory [p^ . The meaning of the pa- 
rameter r can be understood by considering the mean 
{t') = t, and the variance {t''^) — {t'Y = Tt. Hence, r 
rules the strength of time fluctuations, or, otherwise, the 
characteristic correlation time of fluctuations. 

When T ^ 0, Pit.t') 5{t - t') so that p{t) = p(t) 
and V{t) ~ cxp{—iCt} is the usual evolution. 



It is worth noting that the behavior of the distribu- 
tion ( |l^ strongly depends on the regime we consider. In 
fact, for t <C T we have an exponential behavior, while 
for i r a Gausssian-like shape. The case t — t repre- 
sents the border between these two behaviors. All that 
is illustrated in Fig.(|l|). 
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FIG. 1. Probability distribution (^2|) as function of dimen- 
sionless variable t'/r. The dashed line refers to t/r — 0.1, the 
dotted line to t = r, and the solid line to t/r = 5. 

The phase diffusion aspect of the present approach can 
also be seen in the evolution equation for the averaged 
density matrix 'pit). Indeed, by differentiating with re- 
spect to time Eq. and using (|l^) one gets the following 
master equation for 'p(t) 



Xt) = --\og{l + i£T)p{t)- 

T 



(13) 



Once Cp = [H,p], the evolution operator V{t) de- 
scribes a decay of the off diagonal matrix elements in 
the energy representation, whereas the diagonal matrix 
elements remain constants, i.e. the energy is still a con- 
stant of motion. In fact, in the energy eigenbasis, Eqs.(H) 
and (|l^) yield 

Pn,m{t) = CXP (-K„,mi) CXp {-il'n,7nt) /5„,,«(0) , (14) 

where 



Kn,m = ^ log (l + L^I^^t'^) , 
1 

Vn.ni = - arctan (cj„ ^r) , 

T 



(15) 
(16) 



with uJn,m the energy difference. One can recognize in 
Eq.(|l^, beside the exponential decay, a frequency shift 
of every oscillating term. This can be also used as a 
model for Quantum Nondemolition Measurement [p^ . 
In fact, in standard quantum measurement theory each 
measurement results in an instantaneous reduction of 
the wave function onto an eigenstate corresponding to 
the particular detected eigenvalue of the observable be- 
ing measured. Non-selective measurements destroy the 
phase relation between different eigenstates and reduce 
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the state of the system to a statistical mixture where the 
non-diagonal elements of the corresponding density ma- 
trix vanish. Therefore, all random dephasing events, i.e 
all processes that provide for a rapid quantum mechani- 
cal phase destruction but leave the diagonal elements of 
the density matrix unchanged, cause the same dynam- 
ical effect on the evolution of the system like genuine 
quantum- nondemolition measurements. 

However, it would also be possible to consider the Li- 
ouvillian (||) in Eq.jT^), and therefore the competition 
between two types of decoherence. This is what we are 
going to study in the following. 




FIG. 2. Population inversion as function of the scaled (di- 
mensionless) time 7t, for different values of 7r. From bottom 
to top 7r = 0, 0.5, 5, 50. The system has been considered 
initially in the upper level. 



IV. SYSTEM DYNAMICS WITH RANDOM TIME 
EVOLUTION 



If r is small enough, one can expand the logarithm in 
13| ) up to second order in r, and by using the Liouvillian 
, we obtain 



p{t) = -iuj [cr^,p(t)] 
+ 1 {2<j-p{tW 



-o^^ [<^z, Wz,p{t)]] , 



(17) 



where we have used uj ^ j and r ^ 7~^, that is, the 
dissipation takes place on a time scale much larger than 
the time fluctuations. 

Eq.(p^ practically coincides with Eq.(||) provided to 
identify tuP' jl with k. Nonetheless, the present approach 
is different from the usual master equation approach, in 
the sense that it is model independent and without spe- 
cific statistical assumptions. 

For a generic value of r, it is not possible to extract 
an explicit form of master equation from Eq.(^^. Never- 
theless, the physics of the system can be understood by 
simply averaging the quantities of interest over the dis- 
tribution (12). For instance, from Eqs.(p^ and (|^), we 
get 



Tr{p(t)(7j =Tr{75(0)f7jexp 
t 



log (1 7t) 



exp 



log(l + 7t) 



(18) 



Equation (|lg) in the limit 7T ^ 1 reduces to the usual 
decay described by Eq.(||). More generally, the decay 
rate results modified. In particular, for 7r 3> 1 it would 
be possible to inhibit the dissipative effects through the 
nondissipative ones. The frozen dynamics due to increas- 
ing values of 7r is shown in Fig. (|^) . This situation comes 
out as consequence of the transition from the Gaussian to 
the exponential behavior of the probability distribution 
(H) (see Fig.|). 



The freezing effect on the system dynamics remind us 
the quantum Zeno effect Its usual description rests 
on the suppression of the unitary Hamiltonian evolution 
of a quantum system due to intermittent measurements 
in rapid succession. Due to the wave function collapses, 
in the limit of continuos measurements, the evolution is 
completely inhibited, and the system is frozen in its ini- 
tial state. Here, instead, the effect entirely arise from 
quantum statistical properties. 

While in the usual quantum Zeno effect the essen- 
tial requirement is that the measurements of the sys- 
tem state, which cause the interruption, be more closely 
spaced in time than the reservoir correlation (memory) 
time, in our case the correlation time of fluctuations 
should exceed the typical decay time. 

Essentially, we may claim that in the limit 7T 3> 1 
one type of decoherence prevents the other. In fact, 
we may think that dissipative decay process takes place 
through the energy channels determined by the system- 
environment interaction. However, the time evolution 
fluctuations make these channels completely fuzzy, thus 
preventing the decay. 

The above results can be easily extended to the case 
of driven two-level system. 



V. CONCLUSIONS 

In conclusion, we have presented a simple model able 
to explain different aspects of decoherence in a two-level 
system. The used formalism predicts, for large time fluc- 
tuations, a novel Zeno-like effect without invoking the 
abstruse concept of wavefunction collapse |Q . 

The generality of the presented approach suggests in 
some way the possibility that the parameter r (even 
though system-dependent) might have a lower nonzero 
limit, related e.g. to the time energy uncertainty rela- 
tion 1^, or to the finite extension of the spatial wave- 
function ijl^, or even to gravitational effects How- 
ever, even if such "intrinsic" decoherence effects emerge, 
the value of r would be very small. Nevertheless, one 
can think as well to introduce the above statistical prop- 
erties by hand from the outside. For instance, one can 
use dephasing processes through a noisy driving field as 
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envisaged in Ref. ||17|]. Otherwise, one could think at a 
sequence of measurements as jump-hke processes, ran- 
domly distributed in time Q . In such a cases the statis- 
tics, hence the parameter r, would be controlled by the 
experimenter. Thus, it would be an interesting challenge 
to arrange an experimental set up where the conditions 
for above Zcno-like effect are achieved. 
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